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Abstract 

The process p + d <-> 3 He + 7* at intermediate energies is described using a 
covariant and gauge-invariant model, and a realistic pd 3 He vertex. Both pho- 
todisintegration of 3 He and proton-deuteron capture with production of e + e~ 
pairs are studied, and results for cross sections and response functions are 
presented. The effect of time-like form factors on the dilepton cross sections 
is investigated as well. 
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I. INTRODUCTION 



The electromagnetic probe is a well-established and powerful tool to investigate the 
structure of hadronic systems. In exclusive processes one can distinguish three different 
regimes depending on the four- momentum of the photon. First, in the space- like region 
(q 2 < 0) in quasi-free kinematics the (e, e'p) reaction directly probes the single-particle 
structure of the nuclear vertex. Second, for real photons, the (7,2?) or (p, 7) reactions are 
more sensitive to the details of the reaction mechanism and to meson exchange currents 
(MEC). Third, the rather less well-known time-like region (q 2 > 0), which can be explored 
using dilepton production, addresses additional aspects as compared to real photons: (i) the 
coupling of longitudinally polarized photons, (ii) the time-like form factor in the "unphysical" 
region {Am 2 < q 2 < 4m 2 ). Basically there are two alternatives to explore the time-like region: 
virtual Compton scattering p(7, 7*)p and the bremsstrahlung processes with virtual photons, 
such as p + p — > p + p + 7* and capture reactions. 

In this paper we study the capture reaction p + d — > 3 He + 7* at intermediate energies 
(proton energies up to a few hundred MeV) for both real and virtual photons. The main 
motivation for the present work are experiments at TSL (Uppsala) which in the past have 
explored [|lj] only small photon invariant masses (q 2 <(10 MeV) 2 ) and KVI (Groningen) 
experiments which cover larger photon invariant masses 0. 

Since we feel it is important to satisfy gauge invariance, we start from the covariant 
impulse appoximation model of || and make it explicitly gauge invariant by the introduction 
of an additional internal amplitude. An important input in this approach is the pd 3 He vertex 
function, for which recent calculations [ff-|6| using a realistic NN interaction (Argonne vl4 
and vl8) are used. 

II. DESCRIPTION OF THE MODEL 

At small photon energies the photon production amplitude is dominated by radiation 
from the external legs (first three diagrams of Fig. 1), and this consideration led to the devel- 
opment of low-energy theorems (LET)s [0-0. In kinematical conditions where the photon 
energy is not small we may still assume the dominance of the above external amplitude, 
without applying an expansion in powers of the photon energy. 

The 4-momenta of the proton, deuteron, 3 He nucleus, and (virtual) photon are denoted 
by Pi, P2, P3, and q respectively. The external invariant amplitude reads as 

M ext = e e* (A 7 ) u(p 3 , X h ) M^t u(pi, X p ) £ a (A d ) , (1) 

with 

M ext = ® a (P3,P2,Pi - q)S(p! - q,m 1 )T f "(p 1 - q,px) 

+^(P3,P2 - q,Pi)Ap p { P 2 - q)T pa>M (p2 ~ q,p 2 ) 

+r / > 3 ,p 3 + q)S(p 3 + g,m 3 )$ a (p 3 + q,p 2 ,Pi) , (2) 

where e* (A 7 ) and £ a (Ad) are the polarization vectors of the photon and the deuteron respec- 
tively, u(ps, Xh) ( u(p\, X p ) ) is the spinor for 3 He (proton) , e is the proton charge, S(k,m) 
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is the free propagator of the fermion with mass m, Ap p (k) is the deuteron propagator. He- 
licities of the particles are denoted by A's. The structure of the pd— s> 3 He vertex function $ a 
will be discussed later. 

We stress here that in this formulation the s-amplitude (third diagram in Fig. 1) takes 
into account the pole contribution (but not the regular contribution) of the initial state pd 
interaction, and hence takes care of the problem with the orthogonality between initial and 
final state, mentioned in []10[ . 

The electromagnetic (em) vertex function for spin-i particles is chosen in the form 

T^p -q,p) = T»(p, v + q ) = ZY- i—^F 2 (q 2 ) + ZF x {q 2 \q% - gV) , (3) 

2m 

where F\{q 2 ) and F 2 (q 2 ) are respectively the Dirac and Pauli em form factor (FF), F\{q 2 ) = 
[1 — Fi(q 2 )]/q 2 , and Z = 1 (2) for the proton ( 3 He). This vertex obeys the Ward-Takahashi 
identity for the half-off-shell case [|ll| . 



The half-off-shell ^dd vertex satisfying the corresponding Ward-Takahashi identity 
can be written as 



r P°M(p 2 _ QjP2 ) = -g<**(2p 2 - q)» + (p 2 - qfg^ + F 1 (q 2 )g' }a [q 2 (2p 2 - qf - q ■ (2p 2 - q)q»] 
+F 2 (q 2 )(q"g^ - q a 9 w ) + ^^-[q p q a {2p 2 - q)» - \q ■ (2p 2 - q)(q a g^ + g^) ] , (4) 



where F^q 2 ) are related to the charge Gc(q 2 ), magnetic GM^q 2 ) , and quadrupole Gg(q 2 ) 
em FFs of the deuteron (see, e.g., [0). 

Apart from the amplitude M ext corresponding to radiation from the external legs there 
are other more complicated processes, such as initial-state pd rescattering and MEC. This 
contribution (henceforth called the internal amplitude Mi nt ) can be constrained by imposing 
the gauge invariance requirement for the total amplitude M = M ext +M int . Such a procedure 
is conventionally applied in derivations of the LET for bremsstrahlung 0. We will make 
use of consequences of gauge invariance in situations where the photon energy is not small. 

One can show that the internal amplitude obeys the following condition 

q,MZ = -q,MZ = $ Q (p 3 ,p 2 , Pl - q) + $ a (p 3 , P 2 - q,Pi) - 2$> 3 + q,P2,Pi) ■ (5) 

The pd 3 He vertex function for the case where at most one particle is off its mass shell 
has the following structure 

$ Q (fc 3 , k 2 , h) = r + (h, h, h) + r-(h, k 2 , fci A~ mi + ^ 3 ~ m3 C(fc 3 , h, fei) , (6) 

2vri\ 2tti 3 

where the last two terms correspond to negative energy states and only contribute when the 
proton or the helion are off their mass shells. We will use a form for the which allows a 
direct relation with the nonrelativistic wave function (WF) 

r±(h, h, h) = h a G ± (Q 2 ) - Q a H±(Q 2 ) ] 75 , (7) 

where Q a = — ^k 2 is the relative pd 4-momentum and M r = m\m 2 j {m\ + m 2 ) is 

the reduced mass of the pd system. For the 3 He, proton and deuteron diagram in Fig. 1 the 
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relative momenta take the values Q% = ^pf~^P 2 , Qi = Qt~^<f and Q° = Q% + ^q a 
respectively. 

From Eqs. (|5|-[7|) it follows that a solution for the internal amplitude can be constructed 
as MZ = MZ(1) + M&(2), where 

MZil) = { h a G4Ql) - QiH-{Ql)]^- - — [y a G-(Ql) - Q%H-(Q%)] } 7 5 

ZIYI, i 

+g" a [—H + (Ql)-—H + (Ql)h^ (8) 

mi m 2 

M&{2) = —(q - 2pO"J2? + ^(q - 2p 2 )^ 2 Q + (g + 2p S ,)"(J2? + B%) . (9) 

mi m 2 mi + m 2 

We have used the notations ( for 1=1,2 ) 

R? = h«G> + m - Q a 3 H' + (Q 2 )} l5 , G' + (Q 2 ) = g±M^+M) (io) 

and similar notations for H' + (Qf). Note, that this amplitude remains finite in the special 
cases where Q\ — > Q\ or Q\ — > Q\. As a check of our results we verified that M = M ext +M int 
reproduces the LET amplitude of || for real photons when q — > 0. 

We include in the calculation the dominant components [|]|| of the 3 He WF, i.e. a pn pair 
in the deuteron state or in the 1 5'o (quasi) bound d* state, coupled to a proton. We neglect 
contributions to the amplitude where the deuteron is excited into the T = continuum. 
The pd — » pd* capture mechanism via the spin-flip 3 Si + 3 D\ Sq transition has been 
shown || to be important and is therefore included explicitly in M ext . The corresponding 
amplitude (Fig. 1, last graph) can be written as 

MZ = *(p 3 ,p 2 - q, P i)A(p 2 - q)T^(p 2 - q,p 2 ) , (11) 

where A(k) = (k 2 — m* 2 + i0) _1 and the em vertex has the form 

r^(p 2 - q,p 2 ) = -^ u e^q p {p 2 ) v F{q 2 ) . (12) 

Here fi v = Hp — ^ n is the isovector magnetic moment of the nucleon, m 2 is the mass of the 
d*, F(q 2 ) is the transition FF and ^{p^ 1 p 2 — q,pi) is the pd* 3 He vertex function. This con- 
tribution is gauge invariant and does not affect the above discussion of the gauge invariance. 

The invariant functions G±(Q 2 ) , H±(Q 2 ) can be related to the S and D components 
of the overlap integral < d \ 3 He > . For this purpose the formalism developed previously 
in Jl3| for the pnd vertex and also in j|] has been applied. In the same way the vertex 
^(pz,P2 — q,Pi) has been expressed through the overlap integral < d* \ 3 He >. Two models 
for the 3 He WF have been used in calculations. The first one is the parametrization in || 
of the calculations in Ref. with the Argonne vl4 NN + Urbana VII 3N interaction. The 
second model is a more recent calculation |J with the Argonne vl8 NN + Urbana IX 3N 
interaction. 

In order to calculate the em FFs of the proton we used the extended vector meson 



dominance model [14]. The deuteron FFs at negative q 2 are taken from the calculation in 
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||, and we used the parametrization of the FFs of the 3 He as given in |15|]. The proton 
FFs are continuous functions when going from negative to positive q 2 , and this behaviour 
is incorporated in the VMD models. For the deuteron and the helion we have made the 
assumption that the em FFs have a smooth extrapolation from the space-like to the time- 
like region. In this paper we are interested in the interval of relatively small photon invariant 
masses, restricted by the proton incoming energy T p of about 300 MeV. Since the maximal 
photon invariant mass m™ ax = ^ — m 2/3 T p , q 2 does not exceed 0.04 GeV 2 and the 
above approximation should give a reasonable estimate of the effect of FFs on the dilepton 
cross sections. Finally, the FF of the transition d — > d* at positive q 2 is chosen the same as 
the deuteron FF F 1 (q 2 ). 

At the real photon point the FFs are normalized to 

Ff (0) = Tf (0) = G C {0) = 1 , Fi(0) = ^ - 1 , F 2 h (0) = ^ h -2 t 
G A ,(0) = ^/i d , G Q (0) = m 2 2 Q dj (13) 

where Qd = 0.2859 fm 2 is the quadrupole moment of the deuteron and the magnetic moments 
of the proton, deuteron, and 3 He (in nuclear magnetons) are \i v = 2.7928, fid = 0.85774, and 
Hh = —2.12755 respectively. 



III. CROSS SECTION AND RESPONSE FUNCTIONS FOR DILEPTON 

PRODUCTION 

The cm. cross section for the p+d — > 3 He+e + + e~ reaction can be decomposed (in 
complete analogy to the spacelike (e, e'p) reaction) into the sum of products of kinematical 
factors and four response functions (RFs), 

da(e + e-) a 2 mim 3 g c /3 TJ7 1 2 . 2 2 2 



[W T (1- ~P sin 2 9*) + W L (1- (3 2 cos 2 9* 



dQ^dm^dQl 16ir 3 m^p c s 2 

+W TT -P 2 sin 2 9* cos 20* + W LT -^=f3 2 sin 29* cos 0* ] . (14) 
2 2y 2 

Here m 7 = ^fq 1 is the invariant mass of the virtual photon, s = (mi + 7B2) 2 + 2m,2T v , T p is 
the proton kinetic energy in the lab frame, p c and q c are the cm. 3-momenta of the proton 
(deuteron) and photon (helion) respectively, a = 1/137.035, @ — (1 — 4m 2 /m 2 ) 1 / 2 and m e is 
the electron mass. For a description of the kinematics and details about the decomposition 
of the e + e~ cross section into the independent RFs Wj we refer to [ 10| . The differential 



dVt* e in Eq. (|T^) is written in the photon rest frame (denoted by *) [10|, but can easily be 
transformed back to the proton-deuteron cm. frame ( see [|16] ). 

The RFs in Eq. ([14] ) contain information on the hadronic transition. They depend on 
three variables (Wi = Wi(s,m 7 , 9 y )) and are defined as 

1 1 TO 2 

w T = - y2(\j x \ 2 +\j y \ 2 ), w L = -—i v \j z \ 2 , 

6 V 6 Q 2 , 

polar. ™ polar. 

W TT =\Y,{\Jy?-W), W LT = -^ £ 2V2K (J z r x ), (15) 

polar. ^0 polar. 
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where q = (m 2 + q^) 1 ^ 2 is the energy of the virtual photon, and the space components of 
the em current J M = u(p 3 , Xh) M M ° u(pi, X p ) Ca{Xd) are evaluated in the system with the OZ 
axis along the photon momentum. In obtaining these expressions gauge invariance has been 
used to eliminate the time component of the current. 

Integration of Eq. (|14]) over the lepton angles leads to the expression 

da(e+e-) a 2 m im3 q c p(l - \fi 2 ) 

7TV^7 = A^m lPc s [ WT ^ ^' ^ + WL ^ ^' ^ ] ' (16) 

where the interference RFs have dropped out. 

An interesting observable is the ratio of the e + e~ cross section (integrated over the 
allowed photon invariant masses) to the real photon cross section calculated at the same 
incoming energy and scattering angle. This quantity is called the conversion factor, and as 
a ratio it is believed to be less sensitive to many aspects of the reaction mechanism. We 
can cast this ratio in the form R(s, 7 ) = Rt(s,9 7 ) + Rl(s,9 7 ), where the transverse or 
longitudinal conversion factor is given by 

„ ™ max 1 7 

fy r m -y l rim 

RtA*A)= twf n0 J (3(l--(3 2 )q c W T , L (s,m„e,)^. (17) 
irqi Wt(s, 0, tU J2m e 3 m 7 



Here q c = {[s — (m 3 + m 7 ) 2 ][s — (m 3 — m^) 2 ]} 1 / 2 /2^fs , q' c = q c \ mi =Q stands for the real 
photon cm. momentum and m™ ax = ^/s — m 3 . 



IV. RESULTS OF CALCULATIONS AND DISCUSSION 

The model is first tested for the real photon reaction. Fig. 2 (upper panel) shows cross 
sections for the reaction 7 3 He^pd at E lLAB = 245 MeV, related to the capture process at 
T PLAB = 358 MeV via time reversal. Note the unsatisfactory discrepancy between the disin- 
tegration and capture reaction measurements, known for a long time and recently pointed 
out again in Ref . |20| . As seen from the figure the agreement with the photodisintegration 



data [H],nj is quite reasonable and considerably better than with the data [|19|] for the cap- 
ture reaction. The cross section here is determined mainly by the D component of the WF 
and differences between models 'a' ( WF from ) and 'b' ( H ) show up mainly at 9 < 90°. 

To study the importance of the different contributions to the amplitude we plot in Fig. 2 
(lower panel) the energy dependence of the 7 3 He— *-pd cross section at fixed angle. The figure 
shows in particular the large contribution of the internal amplitude (marked 'hit' in Fig. 2 ). 
Also note that the cross section is a result of interference between all contribution, though 
the effect of the 1 Sq is relatively small. All calculations were performed in the Coulomb 
gauge (of course only the amplitude including the internal contribution is independent of 
the photon gauge). The full calculation (solid line) deviates from the data at E 1 between 
100 and 300 MeV, though the disagreement is not too large. We do not show results for 
model 'b'. In general, it gives higher ( up to 40% ) cross sections at small energies, while 
above 50 MeV the situation is reversed. 

We now discuss the e + e~ production in pd capture and present calculations at the proton 
lab energy 190 MeV (corresponding to the kinematics at KVI, Groningen). The invariant 
mass and angular dependences of the response functions are shown in Fig. 3. The choice of 



6 



a forward angle in Fig. 3 (left panel) was made because in these conditions the longitudinal 
response is enhanced. Both Wl and Wlt are comparable in magnitude to Wt- Fig. 3 
shows the sensitivity of the RFs to such ingredients of the model as the time-like FFs, the 
"negative energy" components G_, H_ in the pd 3 He vertex, and the 1 Sq contribution. The 
contribution of the components G-, H is small, which can be explained by a cancellation 
between Mi nt and the terms in M ext proportional to 7 M . Backward angles are less favourable 
for studying the longitudinal response, which is almost always heavily suppresed compared 
to Wt- The exception is at large m 7 , close to the kinematical limit. For that reason we do 
not present the m 7 dependence for backward angles. 

Fig. 3 (right panel) demonstrates the angular dependence at fixed m 7 = 65 MeV (which 
is about \m™ ax ). The transverse RF has a dependence similar to that for the real photons. 
The longitudinal RF is quite large at forward angles and diminishes at backward angles. The 
interference RFs Wtt and Wlt also show up in the forward hemisphere. The interference 
RFs vanish at # 7 = 0° and 180° because of the rotational symmetry around the photon 
momentum in this kinematics. 

From the experimental point of view the cross section integrated over the photon invariant 
masses and the conversion factor are of considerable interest. These observables are plotted 
on Fig. 4. 

Results for R(s, 6* 7 ) (Fig. 4 , lower left panel) are almost independent of the model for the 
WF, despite differences between the cross sections for these WFs (Fig. 4 , upper left panel). 
The effect of the em FFs also turns out to be very small in this kinematics, of the order 
of 1%. In general the longitudinal response is more sensitive to the FFs (see also Fig.2); 
however, after integration over m 7 the longitudinal cross section becomes extremely small 
compared to the transverse one (Fig. 4, upper left panel) except at very forward angles. This 
makes the study of the longitudinal response in this integrated observable experimentally 
difficult. On Fig. 4 we also show a model independent estimate for the conversion factor, 
which follows from Eq. ( ]TT| ) by neglecting Wl{s, m 7 , # 7 ) and the invariant mass dependence 
of the transverse response, i.e., assuming that Wt(s, m 7 , 6* 7 ) ~ Wt(s, 0, 6 I 7 ). On average, the 
deviation between the R(s, # 7 ) in our model and the above estimate are of the order of 5%. 

At higher energy T p = 350 MeV ( m™ ax = 229.8 MeV ) one expects a larger influence 
of the FFs. The differences between the 3 He models 'a' and 'b' in the conversion factor are 
negligible and only calculations with the model 'a' are presented. As it is seen, the FFs 
modify primarily the longitudinal cross section, because the transverse one gets its main 
contribution from the low m 7 region, where q 2 dependence of the FFs can be neglected. 
Compared to T p = 190 MeV, the weight of the longitudinal cross section is now enhanced 
and as a result the conversion factor is more sensitive to the FFs, though the effect is still 
not larger than 5-7%. 

Note, that the time-like FFS of the deuteron and 3 He, when extrapolated from q 2 < 
region, rapidly increase with increasing q 2 . At y/q 1 = 230 MeV, for example, they are 
increased by a factor two compared to their values at q 2 =0. This may not be quite realistic, 
since the extrapolation does not take into account the off-shell effects in the em vertices and 
propagators. The general problem of the (half off-shell) time-like form factors for a weakly 
bound composite system like the deuteron or 3 He is interesting, and will be addressed in 
future work. 

In conclusion, a covariant and gauge invariant approach has been developed for the p+d 



7 



<-> 3 He+7* reactions. The agreement with data for the 3 He photodisintegration is reasonable, 
indicating that the approach seems to account for the basic mechanisms of this process over 
a wide range of energies. An important element of the approach is the internal amplitude 
needed to ensure gauge invariance. The contribution of this part of the amplitude is sizeable, 



which goes in line with the observations made in Ref. ||24|| . 

A mechanism which is missing in this approach is the initial- (or final-) state pd rescat- 
tering, although part of it is effectively taken into account by the 3 He diagram and the 
internal contribution. The effects of explicitly including the pd interaction will be studied 
in a forthcoming publication |J25[| . Predictions have been made for the dilepton production 
experiments under way in Uppsala (TSL) and Groningen (KVI). In general, the longitudinal 
Wl and the interference Wlt response functions are more sensitive to the time-like form 
factors. However, this effect can only be seen at forward angles where Wl and Wlt are 
comparable in magnitude to Wt, or at large angles and high photon invariant masses close 
to the kinematical limit. Finally, we have calculated the conversion factor, which proves to 
be an almost model-independent observable. 
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FIGURES 




FIG. 1. Diagrams corresponding to the external amplitude. 
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FIG. 2. Cross section for the 7 3 He — > pd reaction. Upper panel: angular distribution, lower 
panel: energy dependence. Curves on the upper panel are: solid and dashed ones for the model 
'a' ( WF from || ), dotted and dash-dotted ones for the model 'b' ( WF from Q ). Dashed and 
dash-dotted lines are calculations without Sq amplitude Eq. (TTTT) , solid and dotted lines include 



all contributions. The photodisintegration data are taken from [17 and [18|. The TRIUMF points 
are calculated from the capture data [^] at T p = 350 MeV. Calculations shown on the lower panel 
are performed with model 'a'. The different curves correspond to calculations including different 
contributions to the reaction amplitude. Data are from : o frisl] , V 1 21 1 , • [17|, o [22|, x [23|, and 
* & 
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FIG. 3. Response functions for the e + e~ production in pd capture for proton LAB energy 190 
MeV. Left panel: photon invariant mass distribution at # 7 = 10°, right panel: angular distribution at 
fixed m 7 =65 MeV. Calculations are performed in the model 'a'. Solid lines are the full calculations, 
dashed ones are obtained without the 1 S'o contribution. Dotted lines are the results when the em 
FFs are switched off, dash-dotted lines are results without the G_, components in the pd 3 He 
vertex. 



12 




30 60 90 120 150 180 30 60 90 120 150 180 



e y (cm.) [ deg ] d y (cm.) [ deg ] 



FIG. 4. Angular dependence of the cross section integrated over m 7 (upper panel) and the 
conversion factor devided by the fine-structure constant (lower panel). On the left panels the 
transverse (T) and longitudinal (L) cross sections are calculated with different models for 3 He WF. 
The right panels show the effect of the form factors. The model independent estimate for the 
conversion factor (see text) is shown by the dash-double-dotted lines. 
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